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Abstract. The problem of (non)random distribution of points 
on the sphere is investigated. Published procedures for obtain- 
ing preferred direction and preferred plane for points on the 
sphere (in the sky) are discussed. It is shown that the pub- 
lished methods are incorrect, and, as a consequence, the re- 
sults obtained by these methods cannot be considered to be 
significant. 



1. Introduction 

Neslusan (1995, 1996) and other authors (see references in 
Neslusan 1996), e. g., Tyror 1957, Yabushita 1979, Bogart and 
Noerdlinger 1982 have discussed the distribution of perihelia 
of long-period comets in the sky. The aim of this paper is to 
discuss the mathematical treatment of the methods used by 
these authors. 

2. Perihelion Point Preferred Direction 

Perihelion point of a long-period comet in the sky may be char- 
acterized by direction cosines U, rrii, m. By definition of direc- 



tion cosines, the relation l t + 



+ n- 



= 1 holds. 



If we have N conrets (perihelion points in the sky) , we have 
a set of A triples of direction cosines: {l^fLj, {m,i}fLi, {wijili 
(subscript i denotes the i-th comet). 

Neslusan (1995, 1996) (and other authors - see references 
in Neslusan 1996) finds the preferred direction of the points of 
perihelia and the result is presented by Eqs. (l)-(3) in Neslusan 
(1996): 



in Neslusan's (1996) paper, which is a generalization of the 
preceding papers. However, none of the authors present any 
argument why one should use Eqs. (l)-(3) (with w% — 1 for i — 
1 to N, or, Wi t^ u>j in general). They could be correct only 
if {li}fL 1 ,{mi}fL 1 , {riij^x would be independent quantities 
(the sample mean - the expected value). Moreover, it is only 
the author's decision if he calculates the preferred perihelion 
direction for comets discovered from a given hemisphere with 
weight w = 1/2 for comets discovered independently from both 
hemispheres (the same holds for the older papers w% = 1 for 
i = 1 to N). 

The statement made by Neslusan is that Eqs. (l)-(3) de- 
termine direction cosines < I >, < to > and < n > of the 
preferred direction. This statement is used also in Eq. (5) in 
Neslusan (1996), where it is supposed that < I >, < to > and 
< n > form components of a unit vector (and, the result of 
Eq. (5) is used in Eq. (4) in Neslusan 1996). 

The error is evident. Quantities < I >, < m > and < n > 
do not form components of a unit vector! They are not direction 
cosines of any unit vector (of a preferred direction). Moreover, 
there is no suggestion how to calculate errors for the quantities 
presented by Eqs. (l)-(3). 

3. Other Comments on the Preferred Direction 

According to Neslusan (1996) and other authors, the signifi- 
cance of the departure of the perihelion directional distribu- 
tion from random distribution is characterized by the mean 
quadratic distance of the directions of perihelion points from 
their preferred direction (see Eq. (4) in Neslusan 1996): 
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where Wi is the weight of the i~ th comet in the given calcu- 
lation. It is supposed that some weights may differ from unity 
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The statement is: "In the case of the random distribution, 
< D 2 >— 1/3. The departure from randomness is then char- 
acterized by expression 1/3 - < D 2 >." (below Eq. (5) in 
Neslusan 1996). 

Our comment is that < D 2 >= 1/3 (more correctly, if 
Wi — W, i — 1 to N) holds only for the special type of ran- 
dom distribution - for uniform distribution (uniform in spher- 
ical angles). Moreover, the obtained results for 1/3 - < D 2 > 
Neslusan compares with "the standard deviation", which is 
not defined in Neslusan (1996). It is defined only in Neslusan 



(1995) (Eq. (7), p. 39). However, it does not seem to be correct, 
since for equal weights (the standard case used by authors up 
to Neslusan's work) it yields identical (zero). 

No preferred direction exists for uniform distribution. Since 
the observational data yield a preferred direction, it is impor- 
tant to calculate errors of the investigated quantities. 



of errors can present the significance of the calculated quanti- 
ties. 

Tyror's method (Tyror 1957, used also, e. g., in Bogart 
and Noerdlinger 1982) of calculation of the errors for the sig- 
nificance of the determination of the preferred plane can be 
summarized in the following equations: 



4. Preferred Direction 
Function 



The Probability Density 



A, = {p-Xif , <A>= ^TXi/N 



Eqs. (6)-(9) in Neslusan (1996) define the procedure of com- 
paring the real distribution of the data with the theoretical 
X 2 — distribution (there is not correct sign in the exponent in 
Eq. (9) in Neslusan (1996), and, Eq. (9) does not represent the 
distribution function but the density function) . Unfortunately, 
the procedure used in Neslusan (1996) is incorrect. Neslusan 
defines by Eq. (8) the quantity which is the argument of the 
density function defined by Eq. (9). (The motivation in Eq. 
(8) is incorrect, since the random quantities I, m, n are not 
independent - I 2 + m 2 + n 2 — 1 - and their distributions 
are not of type N(0, 1).) However, theoretical density func- 
tion does not contain the number of the observed elements 
(points), in reality! If one would like to use the idea defined by 
Eqs. (6)-(7), then the correct argument in the density function 
should be the corresponding area of the sphere (e. g., in square 
degrees). Thus, Figs. (3)-(5) in Neslusan (1996) are of no use- 
fulness - it only seems that real distribution may be 
approximated by uniform distribution. However, Figs. 8 
in Neslusan (1995) show that in reality the distribution is not 
exactly uniform - the effect of galactic potential partially dis- 
turbes the uniformity (Neslusan 1995: it seems that directions 
of perihelion points of long-period comets avoid galactic plane). 



5. Preferred Direction — One More Comment 

The section 5 in Neslusan (1996) discusses partial compensa- 
tion of asymmetry in comet discoveries. It seems to us not 
wise to make such a procedure, since we can make whatsoever 

- what we want. It is seen from Eqs. (10) -(11) in Neslusan 
(1996). We have already commented the definition of (original 

- initial) weights. However, Eqs. (10)- (11) modify the original 
weights - the original weights are again weighted in order to 
obtain some better weights. What are the better and the best 
weights? Those, which would produce uniform distribution? In 
any case, any kind of weights yields some preferred direction 
(only special cases yield uniform distribution, however, these 
special cases cannot be easily found and justified). The au- 
thor's role is to explain the obtained result for the existence 
of the preferred direction - it is not enough to say that his 
kind of weights yields different result from that obtained by 
the previous authors. 

The correct procedure should have an opposite direction: 
one must start with the assumption of uniform distribution, 
and, making some known theoretical considerations, he comes 
to some theoretical prediction; the theoretical prediction must 
be consistent with the preferred direction yielded by the ob- 
servational data. 



6. Significance of the Results 

As we have already mentioned, it is important to calculate er- 
rors for the calculated quantities. Only the correct calculation 
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where Xi is unit vector of the perihelion position for the 
i—ih comet, p is unit vector normal to the preferred plane. 
Supposing normal distribution with the mean < A > and 
standard deviation cr < \ > , Tyror's method compares the value 
1/3 — < A > with £r<A>. 

The first argument against the method of Tyror is trivial. 
In reality, unit vector p is determined from the observational 
data ("measurements"), and, thus, it is not exact quantity. 
However, its error is not considered in the calculation of ct<a> 
given by Eq. (5). Thus, real error of < A > is larger than that 
given by Tyror's method. 

The second argument against the method of Tyror is that 
although Eq. (5) seems to be equivalent to calculation of ct<a> 
by the procedure of calculations of errors, the procedure of Ty- 
ror is not standard. Tyror's method is equivalent to calculation 
of the error of A, calculated by the standard error-method pro- 
cedure for measurements, from equation (eigen-value problem 
leads to the cubic equation for < A >) 

A 3 - A 2 + {< x 2 >< y 2 > + < x 2 >< z 2 > + 

< y X z > — < xy > — < xz > — < yz > } A + 

< x >< yz > + < y >< xz > + < z >< xy > — 

2 < xy ><yz>< xz > - < x 2 >< y 2 >< z 2 > = , (6) 

where the errors of quantities < x 2 >, < y 2 >, < z 2 >, < xy >, 
< xz >, < yz > are supposed to be known and are calculated 

as, e. g., a <xy> = \J^ =1 (xiVi- < xy >) 2 /N/(N - 1) . 

Thus, defining Xi = X (x 2 , y 2 , z 2 , xiy t , XiZi, yiZi), < A >= A 
(< x 2 >, < y 2 >, < z 2 >, < xy >, < xz >, < yz >), the 
Tyror's method supposes that the relation 
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(7) 



holds, where partial derivatives are calculated on the basis of 
solution of Eq. (6), and, higher orders in Taylor expansion of Ai 



are neglected. Using the right-hand side of Eq. (7), ct<a>, de- 
fined in Eq. (5), yields the result consistent with that obtained 
by the method of Tyror. Eq. (7) , in reality, corresponds to the 
situation when quantities xj, yj, zj, Xiyi, XiZi, yiZi are con- 
sidered to be independently measured (and higher order terms 
are negligible, i. e., measured quantities are, say, normally dis- 
tributed to a certain degree of approximation). Of course, this 
is not true. Thus, the method of Tyror is not correct. 

Finally, we show the third argument that Tyror's method 
is not correct. In reality, as it was already mentioned, Eq. (7) 
is not exact. It would be not exact even in the case when xj, 
yj, zj, Xiyi, XiZi, y t Zi would be independently measured. The 
reason is that higher orders in Taylor expansion are neglected. 
In reality, one would not need to consider an approximation by 
Taylor expansion since Xi could be calculated from equation 
analogous to Eq. (6) - instead of mean values, values of xj, yf, 
zj , xtyi, XiZi, ytZi should be used. This situation corresponds 
to the eigen-value problem 



C'i Pi — Xi Pi 



(8) 



where matrix components are Cfii = xj, Cn2 = Ca\ = 
Xiyi, Ci\z = Ci3i — XiZi, Ci22 — yi, Caz — C«2 = ViZi, 
Ci33 = zj. This is, of course, geometric nonsense, since the one 
point ( Xi, y iy Zi ) ( xj + yf + zj = 1 ) cannot determine 
a plane, which is, according to Eq. (8), characterized by the 
plane's normal (unit) vector p. In accordance with Eq. (8), one 
would consider the relevant value Xi = 0, and, 



8. Conclusion 

We have shown that mathematical methods used in the prob- 
lem of determination of preferred directions and preferred 
planes for long-period comets are incorrect. Thus, correct pro- 
cedures must be elaborated, and, subsequently, applied to the 
observational data. The results obtained by such correct meth- 
ods must be interpreted, then. All these problems are presented 
in the following parts of this series of papers. 
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As an example we mention that the Tyror's result N = 448, 
< A >= 0.276 would yield <7<a> = 0.013 which corresponds to 
the approximate value 0.014 (Tyror's result) obtained on the 
basis of Taylor expansion defined by Eq. (7). 

We have shown that the results presented in Tyror (1957), 
Bogart and Noerdlinger (1982) are of no significance since the 
method used in these papers corresponds to construction of a 
plane when only one point is known - this is geometrical non- 
sense since, in reality, three points are needed for constructing 
a plane. 



7. Concentration of Points and Its Mathematical and 
Physical Treatment 

Both problems on perihelion concentration of long-period 
comets - perihelion preferred direction and preferred plane for 
perihelion points - were solved by the authors using unit vec- 
tors of individual perihelion points. In general, this method 
yields different results from those obtained on the basis of the 
correct physical treatment of the problem. Physics requires 
that all problems must be solved using radii vectors of the 
perihelion points. One can calculate, using this way, e. g., the 
preferred direction of perihelion points based on the centre 
of mass of the system. Detail theoretical comparison of both 
methods - unit vectors and radii vectors - will be done in the 
following paper. 
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